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Abstract 



Grand unified theories often predict unification of Yukawa couplings (e.g., 
hi, = h T ), and thus certain relations among fermion masses. The lat- 
ter can distinguish these from models that predict only coupling constant 
unification. The implications of Yukawa couplings of the heavy-family in 
the supersymmetric extension of the standard model (when embedded in a 
GUT) are discussed. In particular, uncertainties associated with m* and nib, 
threshold corrections at the low-scale, and threshold and nonrenormalizable- 
operator corrections associated with a grand-unified sector at the high-scale 
are parametrized and estimated. The implication of these and of the correla- 
tion between rat and the prediction for a s are discussed. Constraints on the 
tan (3 range in such models and an upper bound on the i-quark pole mass are 
given and are shown to be affected by the a s — mt correlation. Constraints on 
the low-scale thresholds are found to be weakened by uncertainties associated 
with the high-scale. 

PACS numbers: 12.10.Dm, 11.30.Pb, 12.15.Ff 
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I. INTRODUCTION 



Recent LEP ]I[ and other precision electroweak data is known || to be consistent with 
coupling constant unification within the minimal supersymmetric standard model (MSSM) 
0, in which the standard model (SM) matter is minimally extended, i.e., the Higgs sector 
contains one pair of Higgs doublets and there is a grand desert (up to small perturbations) 
between the weak (low) and unification (high) scales. Recently, it was further shown |4[] that 
corrections associated with the t-quark and Higgs scalar thresholds, sparticle spectrum (for 
example, see Ref. ||), Yukawa couplings, a possible embedding of the MSSM in a grand 
unified theory (GUT) ||, and nonrenormalizable effects 0, as well as constraints ||[| from 
proton decay non-observation |f[0[j , introduce theoretical uncertainties but do not alter the 
successful unification; e.g., the prediction of a s (Mz) ~ 0.125 ± 0.010 M agrees well with 
the observed value. Such uncertainties depend on seven different effective parameters in 
addition to the t-quark mass and Yukawa coupling. (The ±0.010 is a sum (in quadrature) 
of the different theoretical uncertainties estimated using reasonable ranges for the various 
parameters.) This theoretical uncertainty is sufficiently large that few meaningful constraints 
can be derived from the a s (Mz) prediction by itself. Similar conclusions were reached by 
Barbieri, Hall and Sarid |TT| . 

If, indeed, coupling constant unification is a hint for a supersymmetric (SUSY) GUT, 
then a next step is to study the predicted relationships among fermion masses in such 
theories |12||, in a way that consistently incorporates the different theoretical uncertainties 



listed above. (The nature of the theoretical corrections, and in particular the presence of 
adjoint representations, also distinguishes such models from many string-inspired ones.) Let 
us assume in the following (in addition to the MSSM) that we have (i) Coupling constant 
unification, and (ii) Third-family two- Yukawa unification. That is, at the unification point 
Mq (the point above which all the GUT gauge group supermultiplets are complete) we have 
hb(Mc) = h T (Mc), as is the case [12| in a minimal SU(5) unification, which we will assume 
below for definiteness, and in similar unification schemes^]. h a is the MSSM Yukawa coupling 
of a fermion of type a and Mq ~ 10 16 — 10 17 GeV. 

Assumption (ii) can be incorporated into more ambitious attempts (TB| to explain the 
origin of all fermion masses. Such models, which assume extended high-scale structures 
("textures"), were shown recently fT4T]2ll|] to have successful predictions as well as possible 



implications for neutrino masses. However, limiting our analysis to assumptions (i) and 
(ii), we neglect hereafter the Yukawa couplings of the first two families (where empirically 
HM ~ rather than — « — ; the latter would be implied by extending assumption 

(ii) to the first two families and their negligibly small Yukawa couplings) and also flavor 
mixings. The usual argument goes that some perturbation modifies the couplings or masses 
of the first two families^ without significantly altering (ii). We do not elaborate on any 



1 This holds in models such as SU(5), 5O(10) and Eq for the Yukawa coupling of Higgs fields in 
the fundamental (5, 10, 27) representations. 

2 In the texture models mentioned above, such a mechanism is realized by introducing large Higgs 
representations (e.g., 45 of SU(5) or 126 and 210 of S'O(IO)) and (in most cases) a set of flavor 
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such mechanism. A special case of (ii) is a third-family three- Yukawa unification, i.e., 
h t (M G ) = h b (M G ) = h T (Mc), which is the situation in some SO(10) models involving a 
single complex Higgs 10-plet. We will consider such a possibility as well. 

Let us stress that we do not take (ii) to be independent of (i). The coupling constant 
unification assumption by itself is not enough to significantly constrain the MSSM parameter 
space. Here, we examine whether more can be said when imposing (ii) as an additional 
assumption that can possibly distinguish GUT models from some GUT-like string-inspired 
models (where (ii) is not expected to hold in general). Assumption (ii) was considered 
recently by several groups. Some either (a) carried out a one- loop analysis, (b) assumed 
a low a s (M z ) (e.g., a s (M z ) ~ 0.11, which is lower than the value expected from coupling 
constant unification) as an input, (c) ignored the correlation between m t and the predicted 
value of a s (Mz), and/or (d) allowed the running 6-quark mass, m&, to be as high as 5 GeV 
(which, as we discuss below, is a more appropriate upper bound on the pole mass). More 
recent results of two-loop analyses |5,18,|23[] imply a very constrained parameter space, i.e., 



only a small allowed area in the mF t ° 6 — tan 8 plane, where mf ole is the t-quark pole mass and 
tan 8 is the ratio of the two Higgs doublet expectation values, v. Iv . Therefore, one 

wn 

would hope that linking (i) with (ii) (and considering uncertainties associated with m t and 
m&) will result in some useful constraints on the MSSM parameters, assuming a minimal 
S£/(5)-type unification (for example, see Ref. |J). 

Below, we carry out a careful analysis under the above assumptions and consider various 
theoretical uncertainties in the calculation. We find that requiring (i) and predicting a s (Mz) 
as a function of m\ ole and of tan 8 Q in the range of ~ 0.12 — 0.13 (see Figure 1), constrains 
the tan/3 range allowed by (ii) more severely than suggested by previous analyses. On the 
other hand, various theoretical uncertainties can relax the constraints. We also obtain ~ 215 
GeV for the upper bound on rriF t ole (where a s — m t correlations were taken into account). Some 
information about the low-scale mass parameters can be extracted. However, corrections 
associated with the high-scale contribute significantly to the theoretical uncertainties and 
weaken any constraints. The only spectrum parameter that is strongly constrained is tan 8. 
In agreement with other authors, we find low- (~ 0.6 — 3) and high- (~ 40 — 58) tan 8 
allowed regions (branches). The former saturates the ht infra-red fixed-point p4| line (the 
divergence line). The a s — m t correlation modifies the fixed-point value for ht and diminishes 
the dependence of the allowed tan/3 range on vn? t ole < 215 GeV. Theoretical uncertainties 
(and in particular, those associated with the high-scale) determine the width of each branch 
and, thus, the separation between the two branches. 

The various data (and in particular, the 6-quark mass) and the procedure are reviewed 
in section [0. The constraints on the m? t ole — tan 8 plane and the role of the strong coupling 



are presented and discussed in section |HI| . The different correction terms are described 
and evaluated in greater detail in section [TV[ We summarize our conclusions in section 
|V[ Throughout this work, we keep the philosophy (and where relevant, the notation) we 
introduced previously [§J. 



symmetries. For a different possibility involving nonrenormalizable operators, see Ref. plj. 
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II. INPUT DATA AND PROCEDURE 



At the Z-pole, 



M z = 91.187 ±0.007 GeV, 



(1) 



and using the modified minimal subtraction scheme (MS) |25|] the weak angle and couplings^ 
are [|,|26|,[27| 



s 2 (M z ) = 0.2324 - 1.03 x 10- 7 (mr e (GeV) 2 - (138) 2 ) ± 0.0003 



(2) 



a(M z 



127.9 ±0.1, 



(3) 



a s (M z 



0.12±0.01, 



(4) 



where we displayed explicitly the quadratic dependence of s 2 (M z ) on mf oie , which is decou- 
pled from the 0.0003 uncertainty ||. The U(1)y, SU(2) l , and SU(3) C couplings are given 



by 



3 l-s 2 (M z ) 



_ g2 ( M z) aT1 H - 

ai(mz) 5 a(M z ) 'a 2 (M z ) o(M z ) ' miu a 3 (M z ) a s (M z )- 

For the fermion masses, from electroweak precision dataQ we have for the t-quark 



respectively. 



m 



pole 



138l 2 °±5 GeV 



(5) 



for a Higgs mass in the range 50 — 150 GeV, which is appropriate for the MSSM. The pole 
mass is related to the MS running mass, m t , to leading ordeif] in a s by m t = (1 — |^)m£ oZe 



The r-lepton (MS running) mass |JT| is given at the Z-pole by 



m T (M z 



1.7486 ±0.0006 GeV, 



(6) 



which corresponds to (g) and mP T ole = 1.7771 ± 0.0005 GeV |2|. 

The situation regarding the 6-quark mass is more complicated. There are ambiguities in 



the extraction of the MS running mass m^. Gasser and Leutwyler 133] point out that there 



3 A predicted a s (M z ) slightly above 0.13, as it is for a heavy enough i-quark (see Figure 1), 
does not contradict (|j); the a s (M z ) prediction still has a fairly large theoretical uncertainty of 
~ ±0.008. 

4 Slightly more recent data yields |§ m p t ole = 134±|| ± 5 GeV (for m h o ~ 60 - 150 GeV and 
including two-loop aa s m 2 corrections) and s 2 (M z ) = 0.2326 ± 0.0006 (m t free). For our present 
purposes the difference with (|2[) and (||) is negligible. The new data will be incorporated in future 
analyses S. 



5 The next-to-leading correction |3C,31] is ~ 2% (depending on a s ). The leading correction given 

pc 

b 



here is ~ 5%. (See also the discussion of mb/m^ ole below.) We neglect the former (together with 



other subleading mt and mP° le effects) while keeping the latter. 
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is no universal prescription for the relevant scale where a s is to be evaluated, which suggests 
that the extraction of mj is to be carried out case by case, or alternatively, for a range of a s . 
(We will adopt the latter.) Gasser and Leutwyler identify (to leading order in a s ) the running 
mass mft(mfe) with the euclidean mass parameter. This point was emphasized by Narison, 
who offers an alternative definition of m^m^ ' 6 ) pi |. The different definitions introduce a 



scale ambiguity. Another theoretical difficulty may arise from the role of nonperturbative 
effects in the interpretation of potential models^ p3| . The next-to-leading correction to the 



ratio of the MS running mass to the pole mass was given more recently by Gray et al. |p0 
i.e., 

m b = mf e (l-^-12.4(^) 2 ). (7) 

The above comments call for some caution, especially given our aim of exploring the strongly 
constrained m? t ole — tan/3 plane. Let us then adopt a conservative attitude, i.e., 

m b (5 GeV) < 4.45 GeV, (8) 

which corresponds, for example, to m p h ole < 5 GeV, a s > 0.17, and using ([F[). The next-to- 
leading correction term in (|7|) reduces nib, and m? h ole « 4.5 GeV, a s > 0.2 implies m b < 4 
GeV. For example, m? h ole « 4.5 GeV, a s ~ 0.25 gives mj, « 4.0 GeV when neglecting the 
next-to-leading term, and m b J=a 3.7 GeV when using (^). The mj prediction, on the other 
hand, lies (in general) above 4 GeV. Given the above, we do not specify a lower bound 
equivalent to (|j). Also, requiring m{,(4.45 GeV) < 4.45 GeV (which will correspond to 
mb{ m b) — 4.25 ± 0.20 GeV [Q, where we have doubled the uncertainty) is somewhat more 
constraining (e.g., m b (4.45 GeV) — m b (5 GeV) ~ 0.05 — 0.15 GeV - depending on ot s ). 
We use a(M z ), s 2 (M z ), and the r-lepton and t-quark Yukawa couplings, 

h T (M z ) = mAMz) (9) 
rV ; 174 GeV x cos/5 y J 



and 



pole 



MMz ) = ^--^j 174Ge ; xsin/3 . no) 

to predict[] a s (Mz) and hb(Mz), for a definite point in the mF t ole — tan/3 plane. One should 
note that h t depends on m? t ole also via the a s {Mz) correction in ( |10D (and via the «3 con- 
tribution to the running - see below). As we pointed out, s 2 {Mz) depends quadratically 
on mF t ole . Therefore, we neglect all subleading logarithmic dependencies on mP t ole (for a 
discussion, see Ref. ||), including small corrections to (|TDD . We further neglect the error 
bars in ([!]) - (|) and in (||). Also, are all converted to the DR scheme, using the proper 



6 The constituent mass parameter in these models is identified with the pole mass. 
7 i.e., case (b) in the notation of Ref. 
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step functions [|3~5|| . Using two- loop renormalization group equations (RGE's) itera- 



tively, we are able to predict a a {Mz) and h b {Mz) as functions of m pole and tan/?. We take 
100 < mf e < 200 GeV reasonable conservative range, and constrain tan (3 only by 
requiring the Yukawa couplings to stay perturbative, i.e., h a (p) < 3 where M z < // < M G 
and a = t, b, r. (This range can be also justified by requiring two-loop contributions to 
the RGE's to be less than a quarter of the one- loop ones ||18|| .) We then run down using 
three-loop QCD and two-loop QED RGE's to find m°(5 GeV), where the m" prediction 



is that of m b , but without (theoretical) corrections to the RG calculation. 

In any realization of the MSSM, there are small perturbations (order of magnitude of 
two-loop terms) to the grand desert and unification assumptions, as described above. Thus, 
in general, m b = p m" where p ^ 1 is a correction parameter which incorporates the 
uncertainties in the running from Mz up to Mq- Let us stress that in our formulation 
one does not change the MSSM /^-functions to those of the SM at m t or at some other 
effective scale. Rather, leading m pole effects are accounted for in (0), and all other such 
effects determine p \ A point in the m pole — tan j3 is excluded if either h t > 3 (tan /3 < 1 — 2 
and/or m p t ole > 215 GeV), h b (~ h T ) > 3 (tan/3 > 58), or 

m b (5 GeV) = p~ 1 m° b (5 GeV) > 4.45 GeV. (11) 

Incorporating uncertainties associated with m p t ole and Yukawa couplings (in addition to 
the DR conversion step functions) in the numerical procedure^ we have to further consider 
uncertainties associated with the sparticle and Higgs thresholds, high-scale thresholds, and 
Planck-scale nonrenormalizable operators. For simplicity, we will assume that we have 
one heavy (M H 3> M z ) Higgs doublet that decouples with the sparticles, and another light 
(rrih ~ Mz) SM-like doublet that is responsible for all fermion masses^]. We are able to obtain 
an (approximate) analytic expression for p by expanding one-loop expressions around their 
unperturbed values. This will be carried out in section [TV], where we study the different 
contributions to p in GUT models, and estimate p in the minimal SU(5) model. High- 
scale corrections to the coupling constant unification (and not the details of the sparticle 
spectrum) constitute the larger uncertainty. We take 

p' 1 = 1.00 ±0.15, (12) 

which is a conservative estimate derived for reasonable ranges of the various correction 
parameters. Using (O), the exclusion condition ([TT|) reads 



m b {5 GeV) > 0.85mg(5 GeV) > 4.45 GeV. (13) 



8 In the notation of ref. §, Af"™ 6 ^ ", A* op , and ^J ukawa , are all directly incorporated in the 
calculation. 

9 In such a case, SU(2) breaking effects are, in general, negligible above mt. 
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III. THE m p t ole - tan PLANE 



Given the above, we find that assumptions (i) and (ii) allow a low-tan/? branch and a 
high-tan /3 branch. The allowed parameter space is shown in Figure 2, where the narrow strip 
corresponding to three- Yukawa unification is also indicated. The low-tan /3 branch is shown 
in greater detail in Figure 3, where the lines corresponding to p =1 and ht(Ma) = 2 are 
displayed for comparison. The former is, in fact, the h t infra-red fixed-point p4| line, which 
is the /^-divergence line. (h t (Mz) > h{ txed =>- h t (p) ^> 1 for p < M G .) This point was also 
discussed recently in Ref. |T8| , p3l . p ^ 1 only slightly extends the allowed low-tan j3 range. 
It is also interesting to note that constraints from proton decay via dimension-five operators 
would exclude the high-tan (3 branch for p =1 (i.e., tan/5 ^ 4.7 ||). However, once 
correction terms are included, Mq can grow significantly PfllTTI and no useful constraints 
on tan/3 can be derived from proton decay non-observation ||. For comparison, we show 
in Figure 4 the equivalent parameter space with fll"3] ) replaced by 0.85m°(4.45 GeV) > 4.45 
GeV. The allowed tan/3 range is reduced by ~ 0.03 — 0.10 for the low-tan /3 branch and by 
3 — 4 for the high-tan (3 one. Replacing (|13|) with 0.85m°(~ 5.1 GeV) > 4.6 GeV would 
have a similar but opposite effect (i.e., slightly decreasing the separation between the two 
branches). A smaller (larger) uncertainty in (|i~2|) will have an effect similar to the former 
(latter). The p -range estimate (and the m b upper bound) determine the width of each 
branch, and thus the excluded intermediate tan/3 range. Perturbative consistency (i.e., the 
divergence lines discussed above) excludes the very small and the very large tan /3 ranges 
and determines the upper bound on m pole , mf t ole ^ 215 ± 10 GeV, where the ±10 GeV 
uncertainty is due to p top - 

The /^-divergence line eventually becomes approximately parallel to the tan /3-axis (and 
determines the upper bound mf oZe < 215 GeV). Some intermediate values of tan/3 are thus 
allowed for ra\ ole > 200 GeV. (The hb and h t divergence lines meet near the hb = h t line.) 
Otherwise, Yukawa unification at the grand unification scale is ruled outQ if 2.7 ^ tan/3 ^ 
40. Furthermore, the low-tan f3 branch, where ht ~ 1 and which many would consider a 
more natural choice, saturates the fixed-point line and has to be adjusted to a few parts 
in a hundred (a few parts in a thousand, if p 1 = 1) for a given m\ ole (see Figure 3). The 
large-tan (3 branch is more spread and implies, in general, a much lower h t and h t < hb- (h t 
can still be large for a large enough mf°' e , and h t > hb above the three- Yukawa unification 
strip.) While we find no constraints on m\ ole <■ 215 GeV from two- Yukawa unification, 
three- Yukawa unification is ruled out unless 169 ^ m? t ole ^ 196 GeV. (A slightly larger 
range, i.e., m? t ole > 160 GeV, is allowed when one includes corrections to the h t /hb ratio, 
which induce a ~ 5% theoretical uncertainty. We comment more on this point in section |TVj.) 
One expects mutual implications [E29J between the above observations and radiative-breaking 
of SU(2) ® U(l), an attractive feature of the MSSM that prefers h t > h b p|. 

To demonstrate the effect of calculating nib using the predicted a s [Mz) rather than a 
fixed input value, i.e., of associating the Yukawa coupling unification with the rather high 
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For a large tan 8 (v, <C 



some caution may be required regarding the scale at which 



the Higgs potential is minimized and tan/3 is defined, as was pointed out by Bando et al. [22] and 
by Chankowski 
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values of a s {Mz) predicted by a.\ — a 2 unification, we compare Figure 2 with Figures 5-7. 
There, a s (M z ) is fixed (a s (M z ) = 0.11, 0.12, 0.13, in Figure 5, 6, 7, respectively), and thus 
assumption (i) is relaxed; i.e., for a s (M z ) = 0.11 (0.12, 0.13) there is a ~ 7% (~ 3%) split 
between a 3 (M G ) and the a G defined by «i and a 2 . Let us stress that the different corrections 
are not treated on equal footing in this case, because some are included in p while others 
(like NRO's) are absorbed in the fixed value of a s (Mz)- Furthermore, the appropriateness 
of this decomposition depends on which type of uncertainties shift the predicted a s (Mz)- 
(We elaborate more on this point in section |V].) Nevertheless, the comparison illustrates 
that a low a s (M z ) is preferred by m h . The allowed parameter space for a s (M z ) = 0.11 
(Figure 5) is much larger than that for a s {M z ) = 0.13 (Figure 7). For a lower value of a s 
the radiative corrections that reduce hb are diminished, and thus a given h T (Mc) = hb(Mc) 
implies a lower hb(M z ). However, the low value a s (M z ) = 0.11 requires large corrections to 
the coupling constant unification. 

The above discussion also explains the slight differences between our results and those 
of previous analyses. Requiring (i) and using (||) for s 2 (M z ) imply that a s (M z ) grows with 
m\ ole §, e.g., a s (M z ) ~ 0.12 for m p t ole ~ 100 GeV, and a s (M z ) ~ 0.13 for mf e ~ 180 GeV 
(see Figure 1). Indeed, Figure 2 roughly coincides with Figure 6 for the former and with 

pole. , 

Figure 7 for the latter. h t (M z ) in (|T0| ) ~ ^jjjbf> but is diminished by (the vn? t oe dependent) 
a s (M z ). These all affect the balance between positive and negative contributions to the 
Yukawa coupling RGE's (i.e., the fixed-points), and thus modify the hb(Mz) prediction and 
increase the upper bound on mf°' e . 



IV. THE CORRECTION TERMS 

We now turn to a detailed discussion of the correction parameter, p 1 . The coupling 
constant two-loop RGE's are solvable analytically, and it is convenient to write []39f 



— + b i t + 9 i + H i - A; for i = l,2, 3, (14) 



Oli(Mz) «g 

where t = ^Inffj ~ 5.3 is the relevant scale parameter, and olq ~ ^ is the coupling 
constant at the unification point, M G . = 6.6, 1, —3, for i = I, 2, 3, respectively, are 
the one-loop /5-function coefficients; 6>j w 0.7, 1.1, 0.6, for i — 1, 2, 3, are the two-loop 
corrections; Hi are negligible Yukawa coupling two-loop contributions; and the functions Aj 
incorporate all other corrections to the calculation of order of magnitude consistent with 
Q\. In our scheme, ol\ and «2 are inputs. By taking linear combinations we obtain three 
predictions, i.e., 

a s {M z ) = a° s {M z )[a u a 2 , 4 ] + (a°(M z )) 2 A Qfl [A 1; A 2 , A 3 ], (15a) 
— = -Ut*!, a 2 , e u 6 2 ] + A Q [Ai , A 2 ], (15b) 

t = t°[ai, a 2 , 9 X , 6 2 ] + A t [A x , A 2 ], (15c) 
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where we explicitly separated the two-loop predictions with no corrections (Aj = 0) from 
the contribution of the correction functions, Aj. The expressions for A Qs , A Q(3 and A f are 
given in Appendix [A[ 

The integration of the two-loop RGE's for the Yukawa couplings [36| is rather complicated 
and has to be done numerically. To estimate the theoretical correction terms it is useful to 
display the (one- loop) RGE's, i.e., 



dy a 



i=l,2,3 P=t,b,r 



dt', (16) 



where y a = ^ for a = t, b, r; t' = ^ln-^-; and we have omitted higher-order terms. 
bb-i = —-[5, —3, — y, for i — 1, 2, 3, respectively; and bb-p = 1, 6, 1, for (3 — t, b, r. (6 T;i = 
— |, —3, 0; b r .p = 0, 3, 4; b t -i = — 1|, —3, — y; and fr^ = 6, 1, 0.) The balance between 
the negative b^ai and the positive b a -pyp terms determines the infra-red fixed point in the 
Yukawa coupling renormalization flow From (|T6|) we obtain 

JL fa OL (Mr)\^ 

h b (M G ) = h h {M z ) x J] ol)J[ x F b x 6 b x Pb , (17) 
i=l \ a i \ M z) ) 

and similarly for h T (Mo). The af 1, are the one-loop (OL) couplings (i.e., $i = Hi = in 
Substituting instead two- loop (TL) (or input) expressions one has to compensate 
by properly modifying the two-loop correction, Fj, is the correction due to the non- 
negligible Yukawa contribution at one-loop, i.e., Jb a .pypdt'. pb incorporates the theoretical 
uncertainties in the RG calculation. 

From ([17D and the equivalent expression for h T (and assuming hb(Mc) = h T (Ma)) we 
have 

8 10 

m i( M 2 ) = m r (M z) x (sS^Y x f~> x 9- x „- (18) 

where F = Jf, 6 = q 4 -, and p = j-. Setting 6 = p = 1, substituting the one- loop 
expressions for aij and «g, and assuming negligible Yukawa couplings (i.e., F«l) gives an 
exact well-known one-loop expression. F _1 can be estimated analytically for h t ^> hb, h T 
0, i.e., 



(I + II^Mg))-^, (19) 



which gives ~ 0.68 for ht(Mc) ~ 3. In general, however, a numerical analysis is required 
to fully incorporate (F, G) ^ 1. The m°(Mz) that we calculate is given by (|i~8"l) with p = 1 
and numerical values for F and G. 

Before we turn to a rather technical derivation of the correction parameter p , let us 
discuss a simple toy model and point out the ways in which it gets complicated. If the ideal 
desert and unification assumptions hold, then (neglecting two-loop terms) 

* z = 4r + M°, (20a) 
ai (M z ) a° G 
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a 2 (M z ) a 



-o- + M°, (20b) 



— —q- + ht°. (20c) 
a°(M z ) a G 



We use ( p0a| ) and Q20b|) to define a G and i° in terms of the (input) OL\p[Mz)- We now turn 
on the Ai and A2 correction functions and assume that no other corrections contribute to 
p . The coupling constants are now given by 

-J— = -1 + 6 x t° + A QG + 6iA t - Ai, (21a) 
olx[M z ) a G 



-—— = — + b 2 t° + A aG + b 2 A t - A 2 , (21b) 
a 2 {M z ) a G 



1 1 



a s (M z ) a° s (M z ) a ° a. 



-A aa = — + b 3 t° + A aG +b 3 A t . (21c) 



G 



A aa and A t are determined by the condition A aa + biA t — Aj = for 2=1, 2, while (in 
the present approximation) A as is due entirely to the change in a G and t, i.e., — A Qs = 
A QG + b 3 A t . Also, 



1 

o 3 (M G ) ol g 



— + A -g, (22) 



and the a 3 term in (|18D now reads 




l + -[a° s (M z )A as + a G A aG }j. (23) 
We thus obtain (in the toy model) 

p" 1 = e ^[(as(Mz)-aG)^ as -b3a G At}_ ^4) 



In the more general case A 3 7^ and — A aa = A aG + b 3 A t — A 3 , where A 3 = A^ USY + 
^heavy _|_ ^nro ^ QT low-scale threshold, high-scale threshold, and NRO contributions, 
respectively). NRO's (A^ ) modify only the a 3 value and not any RGE coefficients (see 
Ref. |4[] and below) and can be easily incorporated in our toy model, i.e., ([24]) is still correct 
if — A Qs = A aa + & 3 A t — A^ . The high-scale thresholds are more complicated because 
they not only affect a s (M z ) but also change the /3-function coefficient b 3 and the coefficient 
of the RGE for y b (the b b . 3 ) at the various thresholds. The expression for p will be 
derived below. Ignoring for now the threshold changes in bb- 3 the A 3 eavy contribution to p 
no e f K( m z)-qg]a 3 V ^ j e ^ on iy shift in a s (Mz), which affects the entire t' range in (PTB|) , 
is relevant. The effect of the change in b 3 above the threshold is of second order because it 
only affects a small region of the t' integral. Similarly, the leading contribution of A^ USY is 
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p 1 ~ e l;l&s(M z )-a 3 (~TeV)]A§ USY ^ w \y[ c \y j s a second order in small quantities (because it only 
affects a small region of the t' integral) and is therefore negligible. 
Hence, the corrections to gauge couplings lead to 

p- 1 = e f [(a s (M z )-aa)A' as -b 3a GAt]^ ^5) 

where -A' as = — A a< + A% USY = A an + b 3 A t - A^ RO - A 3 eavy includes all the shifts in 



a s (M z ) except those induced by A§ USY . The additional corrections associated with the 
changes in bb-3 at thresholds will be discussed below. 

A different complication is due to the non-negligible role of the Yukawa couplings. F 
is modified when thresholds are decoupled. In particular, once the heavy Higgs doublet is 
decoupled the Yukawa operators and their evolution are modified. (Recall that we assume 
that we have one heavy (M H ^> M z ) Higgs doublet that decouples with the sparticles, and 
another light (m^ ~ Mz) SM-like doublet that is responsible for all fermion masses.) Also, 
h a (Mc) > 1 near either the h t (low-tan/?) or hj, (large-tan f3) fixed points, and the most 
significant high-scale effect of correcting t° — > t° + At is due to the large Yukawa couplings 
and not to the acA t term. We will therefore treat high-scale A t effects (p t ) separately from 
A as and A a<3 effects. A ag will include A 3 contributions which will be partially cancelled by 
decoupling thresholds from both the a 3 and y b RGE's. Thus, A ag and A ac effects will be 
described by p a:j , which we derive first. (Using the input value of ai, p ~ 1 - see below.) 
We will then consider corrections to F {p^)- Lastly, we will derive p t 1 and rewrite p 1 in 
a way that reflects the correlations among p a3 , p 1 and p t 1 . We will also comment on the 
role of the high-scale corrections, the case of using a s (Mz) as an input, and on corrections 
to the h t /h b ratio. 

Allowing a complicated threshold structure near Mz (and/or near Mq) gives a modified 
one-loop expression for rnj,, 

3 n 

m b (M z ) = m T (M z ) x J[ ]J ( I xr'xfl + A,), (26) 

i=l k=0 \ a i\P ) / 

where k runs over the various thresholds, i.e., p° = M z and p n = Mq- b k is the one-loop 
coefficient of the respective RGE between p k and p k+l ; and represents the threshold 
corrections to F. By expanding ([26|) around ( |T5D (in a similar way to (J£3|)) and using the 




results of Ref. we can obtain an approximate expression for p . This yields a better 
insight into the role of the different correction parameters than purely numerical estimates. 

The important effects of the coupling constant uncertainties are in the a 3 terms. a2 (in 
our approximation^) drops out from (26) and the residual uncertainties from cti are small 



when the input value is used. Recall that our strategy is to use the experimental values 



Once sparticles are decoupled the degeneracy among various operators is lifted, e.g., the gaugino 
- sfermion - fermion coupling is different from the respective gauge coupling and the higgsino - 
sfermion - fermion Yukawa coupling is different from the Higgs-boson - fermion - fermion one (see, 
for example, Chankowski |37fl ). In ( |26|) we ignored this effect, which is negligible for sparticles and 
the Higgs-doublet below the TeV scale. 
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of aci(Mz) and a 2 (M^) to predict 03. The dominant corrections to the ra& prediction are 
the uncertainties in etc and t due to Ai and A 2 , and the explicit uncertainties in A 3 (as 
was illustrated by our toy model). The latter can be divided into low-scale (A§ USY ) and 
to high-scale (Aj™" 9 + A^ RO ) contributions. The low-scale uncertainties have only a small 
effect on m& because they only affect a small t' range in ([16]) (see the toy model). High-scale 
corrections affect the entire if range. They modify both a s [Mz) (high-scale contributions 
to A 3 constitute a part of A Q J and either the (3% function near Mq (Ag 611 ^) or the a^Mo) 
value (A^ ). All (high- and low-scale threshold) corrections to (3% affect the a% terms in 

We denote the heavy X and Y vector; color-triplet; and the adjoint color-octet, SU(2)- 
triplet (and singlet) superfield thresholds by My, M 5 , and M 24 , respectively. Some of the 
high-scale thresholds are strongly constrained by proton decay, i.e., in the minimal SU(5) 
model (which we assume) M 5 ~ Mq and perturbative consistency constrains Mq < 3My 
M 24 <C Mq is possible, and A t in this scenario can be ~ +0.5 and the constraints on 
M 5 are relaxed (i.e., M 5 > O.IMq) ||. Also, proton decay constraints can be removed by a 
simple modification of the model |41| . 

The sparticles and the Higgs doublet decouple from the aij RGE at an effective scale, 
Mi, defined in Ref. M (see also Carena et al. P3[), i.e., 



^ (2ir) M z ~ (2n) M z '° r * " ' 6 ' { "> 

The summation is over all relevant thresholds, i.e., sparticles and the heavy Higgs doublet, 
and b\ is the C-paxticle contribution to the respective /3-function. Mj can be split by a 
factor of a few. In general, M\ grows most significantly with the scalar mass; M3 with 
the gaugino mass; and M\ and M 2 grow the same with the higgsino mass; and M 2 M\ 
and/or M 2 M3. Mi, M 2 and M3 all appear in A Qs , A QG and Af. On the other hand, once 
either the gluinos or the squarks are decoupled, all squark - gluino loops are eliminated and 
bb,3 — ||3711 1 and two other scales of relevance are (in the approximation of degenerate 
squark masses) M 3 = min (M gluino , M squark ) and M 3 = max (M g i uino , M squark ). One has 
M| = M3M3. 

We consider high-scale thresholds and NRO's (~ [a s (Mz) — «g]A), low-scale thresholds 
(~ [a s (Mz) — «3(~ TeV)]A), and corrections to the coupling constant unification predictions 
for a s (Mz) and olg- We will discuss corrections to F and to t below. A more detailed 
treatment of low-scale effects will be needed if either some of the spectrum parameters are 
better known or if one assumes sparticle thresholds above the TeV scale. We will takeQ 
/i 1 = M 3 , fi 2 = M 3 , /i 3 = M 24 and = M 5 . The couplings and coefficients (to be 
substituted in (|26|)) read 

a 3 (/i°) = a° s (M z ) + (a° s (M z )) 2 A aa , (28a) 



The generalization to M5 < M 24 is straight forward. The My < Mq case is much more difficult 
to describe. The heavy X and Y supervectors couple to the SU3 x SU2 X Ui Yukawa operators in 
a complicated way. However, My > \Mq and the effects cannot be large. 



12 



M//))- 1 = {oP s {Mz)Y l - b% - A as , (28b) 
(asifi 2 ))- 1 = (a^Mz))- 1 - b% - b\8h ~ A Qs , (28c) 



(asin 3 ))' 1 = M/i 2 ))- 1 " &3 ^ = = M/x 4 ))" 1 (28d) 



M/))" 1 = (c^)- 1 - b% + A aG - A™, (28e) 
M/* 5 ))" 1 = M/Or 1 = (og.)" 1 + A QG - A™, (28f) 
6° = bf M = -7, 6 3 = -5, bl = bf SSM = -3, b\ = 0, b* = 1, (28g) 

= &£s = C = ~ 8 ' & ';3 = & ';3 = = b^ SM = -y • ( 2 8 h ) 
*3 = ^ ln l§> ^3 = ^ ln ^, and i 5 = ^hig|. (We replaced by the two-loop a° G 



which introduces a negligible inconsistency.) In the [M24,Ms] interval we cannot use (p6| 
because 6| = 0, and instead we have & 4 ;3 a3(//) din//, which contributes — |^f[l n ^7 _ 
ln^]=« G [-|Ar + It 5 ]tolnp- 1 



M G J - "«L 9^3 -r 3 l 

We obtain (for the a 3 terms in 




«3 ' 



(29) 



where 



- 1 ^Q°A _20 t 4 0g t 8 0, 8 a A 24_8 aJVHOi . . 

5^ = gl9 a s Q a s t 9 a G a aG g O s I3 9 O 3 0I3 g a G r 5 9 Q! G ii 3 9 Q G ii 3 I (30) 



(a° is a°(M*).) 

«i (and a^) uncertainties feed into A Qa , A aa and A t (we discuss the latter below). 
There are also p corrections from thresholds and NRO's analogous to ( p0[) from the 



(ai(M z ) /ai(M G ))99 factor. However, these are negligible (< 1% or p ~ 1) when the 
experimental input value for a 1 (M z ) is used. We take in (^Bj) /j = 1. The «i term in ([TO]) 
does, however, lead to a small contribution to the p t 1 term. 

A Qs , A QG , A^ RO and A 3 4 are defined in Ref. Q and are given in Appendix |A| for 
completeness. They involve the low- and high-scale mass parameters introduced above, as 
well as the NRO effective strength, rj. To leading order rj is the only NRO free parameter 
and it incorporates the degrees of freedom associated with the strength, sign, scale, and 
normalization of the dimension-five operators — -Tr(F flv ^F IJ ' I/ ), where F^ u is the field 

strength tensor and $ is the adjoint scalar field. The range —10 < r\ ^ 10 suggested in Ref. 
H is constrained only by perturbative consistency of the analysis. 
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Threshold corrections also affect the one- loop contribution from the Yukawa sector, i.e., 
F — > F(l + Ap), and it is convenient to define 

p- F l = l + A F . (31) 

F^ 1 is a correction term, but it can be as large as a ~ 30% correction (which, in fact, is 
responsible for the successful m& prediction in the MSSM), and, as we shall show, Ap m 
2% — 4%. M 2 4 <C Mq will not contribute since the adjoint superfield couples (to one-loop) 
to the Yukawa operators via its coupling to the Higgs doublets, which drops out from the 
ratio. However, new and large Yukawa couplings will (radiatively) increase h a (ji) and thus 
affect the infra-red fixed points and the perturbative limit; i.e., they affect F a rather than 
the ratio F. (Such an effect may shift the ht and hj, divergence lines in Figures 2 — 7 inwards 
towards each other.) New Yukawa operators (that do contribute to the raticQ) are also 
generated if M 5 < M G (see, for example, Hisano et al. ||). The exact magnitude of such 
effects will be determined by the details of the high-scale Lagrangian. 

There are, however, low-scale corrections to F~ x . We naively change the Yukawa cou- 
pling RGE's below the heavy Higgs doublet threshold (i# = In j^) to those which are 
appropriate given the SM fermion spectrum with one SM-like Higgs doublet (for example, 
see Giveon et al. p2|). We will also neglect (near Mz) hb cos (3, h T cos[3 ~ 0. We obtain 

p- 1 _ e [\yt+\yt sir? i3]t H ^ ^32) 

where here y t is taken at Mz (or more correctly, between Mz and Mh), and tn < 0.38. p 1 
increases m;, by slightly diminishing the effect of F~ l in (p!8|) . Note that the F behavior 
distinguishes the MSSM, where only hb gets corrected (to one-loop) by h t , from the SM 
where all fermions couple to only one Higgs doublet, and both hb and h T get corrected. 
For ht < 1.1 (as is reasonable at the low-scale) p p ^ 1.04, which is a naive overestimate. 
Including a hb(Mz) — h T (Mz) (~ 0.4) contribution can increase p 1 by less than ~ 2% (the 
upper bound is for a large tan/3). In most parts of the plane the correction is moderate, 
i.e., p ^$ 1.02 if either sin/3 ~ or /i t < 1. Let us stress that this is a somewhat naive 
description which gets complicated in many ways. For example, a light t-squark and a light 
chargino will still couple to the SM-like effective Yukawa operators. Such effects will have to 
be accounted for if and when the spectrum is better known and a refined analysis is required. 

Lastly, t (which is determined by a% — ot2 unification) can be corrected by either correc- 
tions to the coupling constant unification (see Eq. or by a split between the coupling 
constant and Yukawa coupling unification points. In the latter case, from our definition 
of M G , A t < (and it is reasonable to take A t ^> —1). (Effects (e.g., NRO's) that may 
split hb(Mc) and h r (Mc) can be also expressed in terms of the split between the unification 
points, but then A t has no fixed sign.) Taking the approximation that At t so that 
ati(t) w an(t + A t ), h t {t) mht(t+ A t ) we find 

p" 1 = e (->+2a° G )A t; (33) 



13 Their effect can be estimated by observation of the SU(5) invariant operators, i.e., Fb/F T — ► 1 
above M5, ((l^) is slightly modified for M5 < Mq, and the divergence lines move slightly outwards. 
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where y t in ( |3"3"D is taken at M G , i.e., y t (M G ) = h * ^ G>> , and h b « /i r dropped out. For small 
values of h t (M G ), a longer running time reduces h b (Ma) (and thus, increases the predicted 
hb(Mz)) and vice versa. The situation reverses for h t (Mc) ^ V^16™^; ~ v2- 
Collecting our results, we have 



where 



p" = Pal >< pJ >< Pt 1 = P^ X Po 1 ' ( 34 ) 



, w x [iCi+^C 2 +15C 8 ] M [fC 1+ f§C 2 +25C 8 ] 

P* = TT 




M z 

ImJ vmJ ' 1 ; 



and 



Pg = VMgJ X VMg / 

x (^j x (1+ [0.29^ + 0.24^ + C 8 ]i7), (36) 

represent the low-scale and high-scale corrections, respectively. The coefficients C{ are de- 
fined and estimated in Table |I[ Note that C\ + C3 — C4 = 0, i.e., M3 drops out. This is 
because M3 is associated with the change in a s (M z ) due the threshold, which is a second 
order effect (see the discussion above). The M 3 dependence, on the other hand, is due to 
the change of the 6^3 coefficient and is of first order. We used Mf = M3M3 and added A aG 
and ^-terms. 

It is instructive to rewrite (using Table |) 

-, fM,^'" ( M,y'' K (WA 0025 fM H \° m 

p * " [wj {wj lm) feJ ■ (37) 

(Different values of CV and Cg were averaged.) If the spectrum were all degenerate at 
Msusy, then p^ ~ ( M ff z sy ) ~ 0.94. We can invert the logic and use (|37 ) to define an 
effective scale that gives p 1 correctly. For example, in Ref. |4[] we defined an effective scale 
parameter, A SUSY , 

Mi M 2 Ms A<?r/<?y 

251n— -- lOOln— - + 561n — - = -191n-^. 38 

M z M z M z M z 

(Asusy here is Msusy of Ref. |§ , and we have changed notation in order to avoid confusion 
with other definitions of Msusy-) Asusy gives correctly the corrections to the a s (Mz) (or 
s 2 (Mz)) prediction, but does not contain any information on the spectrum - it can be as 
low as a few GeV for sparticles 3> M Z - (See also Carena et al. p3"|.) Here, we can similarly 
define 

-1 (Bsusy\ [ -^-^ CM+C ?- WCs] (BsusyY 0025 ,, q . 
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The slightly negative exponent implies in many cases (for non-degenerate spectra) Bsusy ~ 
M z (i.e., M 2 < M%, M 3 , Mh usually implies p z ^ 1). Bsusy > Mz for a strongly de- 
generate spectrum. For the spectra of Ref. |J we find p 1 ~ 1 (A susy ~ 32, 21 GeV, 
Bsusy ~ Mz)- Taking the limits of heavy gluinos and of a degenerate spectrum we find 
0.94 < p < 1.06. Away from the limits p —>■ 1. 
Similarly to (|37|) we can rewrite 

/ !\ yf \ -0.030 , Tiir -0.004 , M N 0.015 

A scenario in which M 24 < M G ; M 5 ~ (0.1 - 0.5)M G ; My = M G ; and rj « -10; would give 
p c « 0.9. This scenario is also consistent with limits from proton decay [§,0. Furthermore, 
NRO's contribute only negligibly to A a<3 and to A f (unless one allows NRO effects to be very 
large fl7|f42"ll ) . M 2 4 <C M G on the other hand can increase t significantly, i.e., M G <5x 10 17 
GeV (which is the reason that we can have M 5 < Mq). A large negative rj maintains an 
acceptable value of a s (M z ) in such a scenario. Lifting proton decay constraints (e.g., see 
Ref. |41[]), we can have M 5 <C Mq and p G ~ 0.8 — 0.9. Taking these limits and that of a 
degenerate spectrum and a large positive rj we obtain 0.8 ^ pj ^$1.1 

The high-scale corrections to the coupling constant unification emerge as the leading 
contribution to p ^ 1. We would like to stress that rj is not just a new ad-hoc parameter. 
Given the precision to which we know the low-scale observables, one cannot ignore the likely 
possibility of unknown physics at the high-scale where the (supergravity-induced) MSSM 
breaks down, and which is parameterized in terms of NRO's (whose form is defined in SU(5) 
models). Furthermore, similar corrections may arise in supergravity from non-minimal (and 
non- universal) gauge kinetic functions (see, for example, Ref. [42]). Unfortunately, this, in 
turn, introduces some ambiguity in RG calculations (via high-scale boundary conditions). 
It should also be noted that adding large representations [p!3] - [T6|j20| , e.g., 126 of 5*0(10), 



does not introduce (for nearly degenerate heavy components) large threshold corrections 
to a s (Mz) and t. This is because the decoupled heavy components constitute a nearly 
complete representation (which acts equally on all the Vs). Thus, the threshold corrections 
in the minimal model give a good estimate of p G (in models with a GUT sector, which 
are the relevant ones for Yukawa unification). A model independent treatment of high-scale 
threshold effects on coupling constant unification was given in Ref. 0. The heavy Yukawa 
sectors of different models may affect the infra-red fixed-points differently. 

An arbitrary splitting of the two unification points induces a ~ 5% uncertainty. By 
combining all the contributions in quadrature (as a guideline only) we obtain 

0.80 <p~" < 1.15. (41) 

p = 1 ± 0.1 is thus a reasonable range, and p 1 = 1 ± 0.15 (that we adopted) is somewhat 
more extreme, but well within the allowed range. We would like to stress that all ranges 
extracted here are a guideline only. This range, which is controlled by high-scale correc- 
tions, is still valid when the sparticle spectrum is explicitly calculated (and, e.g., decoupled 
numerically) . 

As we pointed out above, corrections that either change the prediction for a s [Mz) o r the 
positive contribution to (|T6j) from Yukawa terms, affect the infra-red fixed points and can 
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slightly shift the corresponding divergence lines in Figures 2 — 7. They may also affect the 
upper bound on mf°' e , and thus they induce a ~ ±10 GeV uncertainty to the upper bound 
value. However, the corrected range remains > 200 GeV (see also Barger et al. |18[). In 
particular, it is much higher than the the upper bound suggested by precision data (see Eq. 

(!))• 

We would also like to point out that if ai(Mz), a^iMz) and aj,{Mz) are all used as 
inputs, then one arbitrarily adjusts A Qs so that oP s + (a°) 2 A Qs is fixed at some desired 
value. The coupling constants do not unify unless one consistently corrects «g and t as well. 
However, such a procedure is a reasonable approximation if Af f/5Y is small (or is known 
and corrected for). In that case one can minimize the residual uncertainty by calculating 
«3(Mg) from the input value of a s {Mz) and from Mq [Mq is calculated from a\ — a<i 
unification). Then only p t , p p and bb-3 terms contribute to p (i.e., one can obtain their 
contribution by setting G\ = C 2 = in (^) and (p6|)), and the residual uncertainty is small. 
Some caution is, however, needed. The coupling constant unification constraints are not 
integral in such a procedure (e.g., compare Figures 2 and 5). In particular, the correlation 
between a s (Mz) and mt is not manifest. A large m t value implies larger values of a s (Mz) 
or, alternatively, very large corrections to coupling constant unification. Also, only Ai, A2, 
A^ RO and /\^ eavy can induce first order corrections to mj, and thus can be used to fix 
a s (Mz)- (NRO's renormalize and split a^Mc), and thus honestly modify the boundary 
conditions. The simplest way to adjust the a s {Mz) prediction to a given input value is by 
adjusting 77 - i.e., r\ ~ —10 corrects the a s {Mz) prediction to a s {Mz) ~ 0.11.) As was 
illustrated by our toy model, Af USY contributes to A aa but does not affect mj to first order 
in small terms. Thus, unless one knows and corrects for the A§ USY contribution to the 
input a s (Mz) one introduces a significant theoretical uncertainty. Lastly, the experimental 
uncertainty in a s (Mz) is large, and arbitrarily varying a s [Mz) i n that range is not very 
instructive. Nevertheless, it is useful in demonstrating the role of a s (Mz) in predicting rrib, 
as we saw in section |TJ. 

Finally, the three- Yukawa unification strip (see Figure 2) has uncertainties in both the 
tan/3 and the m? t ole ranges, coming from corrections to the hb/h T and h t /h b ratios, respec- 
tively. To one-loop (h t /hb) ~ (af i /a§ i )~33 x (Fb/F t ) and any uncertainties in the ot\ term 
are negligible. However, variation of —0.5 < A t < 0.5 generates ~ ±2% (~ ±8%) correc- 
tion if h t ~ h h ~ h T ~ 1 (~ 2), i.e., (pb/pto P )t ~ e 2^( M e)A t _ Additional uncertainty of 
0.95 < {pb/pto P )F ~ e (j/t(Mz) "t^ (Mz)) * ff < 1 (we assume tan/3 > 1) is associated with the 
decoupling of the heavy Higgs doublet. We estimate a ~ ±5 — 10% uncertainty in the mF t ole 
range that corresponds to three- Yukawa unification. 



V. CONCLUSIONS 

Grand unified theories typically predict h b = h T at Mq, and contain non-fundamental 
Higgs representations. These distinguish such models from some other realizations of the 
MSSM, e.g., string-inspired GUT-like models. Above, we explicitly embedded the MSSM in 
a minimal SU(5) model, and concluded that such a model is constrained to a small area of 
the parameter space. We showed that corrections to a two-loop calculation of the bottom 
mass (when assuming grand unification) are manifested in various ways. Parametrizing 
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those corrections, we were able to relate them to the correction parameters identified in Ref. 
ID, and to study their magnitude and behavior in some detail. The theoretical uncertainty 
in the bottom mass prediction is typically < 15%. We thus took (given the ambiguities 
in the extraction of mi from experiment) 0.85m°(5 GeV) < 4.45 GeV as a (conservative) 
constraint. Requiring this, as well as requiring perturbative Yukawa couplings up to M G 
and identifying the coupling constant and the (third-family) Yukawa coupling unification 
points, we found that the range 2.7 ^ tan/5 < 40 is excluded (as well as mf le ^ 215 GeV), 
and that, in agreement with other authors, the allowed area in the m p t ole — tan/3 plane is 
described by low- and high-tan (3 branches (where the former saturates the h t fixed-point 
line). The separation between the two branches is determined by the correction factor. 
Requiring all three (third-family) Yukawa couplings to meet constrains 160 GeV ^5 mF t ole 
and requires a large tan/5. We demonstrated that the allowed parameter space grows for 
lower (input) values of a s (M z ), but that the MSSM prefers higher values. We further 
argued that the s 2 (M z ) quadratic dependence on wP° le cannot be ignored as it correlates 
the a s (Mz) prediction with mf°' e , and thus affects the mf 3 ' 6 dependence of the prediction, 
as well as the the upper bound on m\ ole and the range of m\ ole for which intermediate values 
of tan f3 are allowed. Finally, we expect the above observations and radiative breaking of 
SU(2)®U(l) to have mutual implications, and suggest that the above constraint is still valid 
in a calculation in which the sparticle spectrum, and therefore p , is calculated explicitly. 
(The larger uncertainty in the calculation comes from the unification- scale physics rather 
than from the details of the sparticle spectrum.) Our hope is that a careful study of various 
correction terms will eventually result in reliable constraints on the MSSM parameter space, 
and in a way that can distinguish different realizations of the MSSM. Here we have showed 
(in agreement with others) that by measuring tan/3 one can exclude simple (and some 
extended) GUT structures at the high-scale. 



ACKNOWLEDGMENTS 

This work was supported by the Department of Energy Grant No. DE-AC02-76-ERO- 
3071. 



APPENDIX A: THE CORRECTION FUNCTIONS 

For completeness, we give the correction functions to the coupling constant unification 
(in the minimal SU(5) MSSM). For more details, see Ref. ||. Corrections that depend on 
rnF t ole or the conversion to the DR scheme are included in the numerical procedure, and are 
not quoted below. All the parameters are defined above (see section [TV]) . 

, My , M 24 , M 5 
28ttA q = -12 In— - - 61n— p + 18m — - 

Mq m g m g 

, M, M 2 , M 3 , k . 

+25 In— - - 100 In + 56 In-/ + 8.00/7. Al 

M z M z M z ' ( J 

- 336vrA ar = +888 In — ^ - 396 In — + 12 In — 

M G M G M G 
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+751n^- - 825 In ^ + 50.0^. (A2) 
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m 12 m 1 In 

Mq M g 5 Mq 

M. , M 2 . . 
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FIGURES 



FIG. 1. The predicted strong coupling at the Z-pole, a s (Mz), for different values of the i-quark 
pole mass, mf°' 6 , and of the two Higgs doublet expectation value ratio, tan/3. h^Mo) = h T (Mc) 
is assumed. mf t ole (in GeV) is indicated on the right-hand-side above the relevant line. 

FIG. 2. The mf°' e — tan/3 plane is divided into five different regions. Two areas (low- and 
high-tan (3 branches) are consistent with perturbative two- Yukawa unification (h^Mo) = h T (Mo)) 
and with 0.85ra^(5 GeV) < 4.45 GeV. Between the two branches the 6-quark mass is too high. For 
a too low (high) tan/3, h t (hi,) diverges. The strip where all three (third-family) Yukawa couplings 
unify intersects the allowed high-tan (3 branch and is indicated as well (dash-dot line) . Corrections 
to the ht/hb ratio induce a ~ ±5% (vertical) uncertainty in the m\ ole range that corresponds to 
each of the points in the three- Yukawa unification strip. a s (Mz), ag, and the unification scale 
used in the calculation are the ones predicted by the MSSM coupling constant unification, and 
are sensitive to the t-quark pole mass, m? t ole (see Figure 1). The wP t ole range suggested by the 
electroweak data is indicated (dashed lines) for comparison. mf ole is in GeV. 

FIG. 3. The low-tan (3 branch of Figure 2 is shown in greater detail. The lines corresponding 
to p =1 (thick) and ht(Mc) = 2 (dashed) are indicated. To the left of the allowed branch one 
obtains ht{Mc) > 3. 

FIG. 4. The same as Figure 2, except the constraint is replaced with the more restrictive one, 
0.85m£(4.45 GeV) < 4.45 GeV. The allowed tan/3 range is reduced by ~ 0.03 - 0.10 for the 
low-tan (3 branch (the effect is hardly seen in the figure) and by ~ 3 — 4 for the high-tan (3 branch 
(where the corresponding range for 0.85m°(5 GeV) < 4.45 GeV is indicated - dashed line - for 
comparison). 

FIG. 5. The area in the wP° le — tan (3 plane which is consistent with perturbative two- Yukawa 
unification and with 0.85771^(5 GeV) < 4.45 GeV assuming a s (Mz) = 0.11. The unification scale 
and «g used in the calculation are those predicted by a\ — Q2 unification. We chose p = 0.85 
for comparison with Figures 2 — 3. mP t ole is in GeV. 

FIG. 6. The same as Figure 5, except a s (Mz) = 0.12. 

FIG. 7. The same as Figure 5, except a s (M z ) = 0.13. 



22 



TABLES 



TABLE I. The coefficients Q are defined and estimated using s 2 (M z ) = 0.2324, a° s = 0.125, 
and a% = 0.040. 
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estimate 


comments 


Ci 


8^>(M Z ) 

9 TT 


+0.035 




c 2 


o Q_ 

9 7T 


-0.011 
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10 oP s {M z ) 
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-0.009 




c 5 
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3 7T 
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M 24 < M 5 
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9 n 
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2+3 sin/3j/t(M z ) 
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+0.010 
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-0.0013 
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